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Abstract 
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I. INTRODUCTION 



Properties of hadrons at high temperatures and densities are of great interest from both 
experimental and theoretical perspectives. One motivation for studying temperature ef- 
fects on hadrons comes from the suggestion that some features of the dilepton spectrum 
observed in heavy-ion collisions can be explained by the modification of masses and widths 
of mesons by the thermal medium ffl]. Nevertheless, reliable information on the temper- 
ature modification of hadronic properties is still lacking. Lattice simulations, which rely 
on the imaginary-time formulation of quantum field theory, have serious difficulties with 
real-time quantities. 1 The absence of Lorentz invariance at finite temperature implies that 
there is no direct relationship between real-time characteristics of hadrons (for example, the 
so-called "pole masses," which are supposedly the positions of poles in propagators) and 
quantities that can be extracted from Euclidean propagators (e.g., the "screening masses," 
which characterize the exponential falloff of static Euclidean correlators). Thus, as a rule, 
lattice measurements of correlation functions at finite temperature cannot be used to draw 
conclusions about real-time propagation of hadrons. 

The aim of this paper is to demonstrate that pions present an exception to this rule. We 
shall argue that it is possible to determine the dispersion relation of soft pions (more precisely, 
its real part) at all temperatures below the chiral phase transition, knowing only equal-time 
(or static) correlation functions, which, in principle, can be determined on the lattice. It 
should be emphasized that we do not assume the temperature T to be small compared to the 
chiral phase transition temperature T c : we must have T < T c , but T/T c is allowed to be of 
order 1. Our results are valid in the nontrivial regime where neither perturbative QCD nor 
chiral perturbation theory are reliable. Moreover, our method, with minimal modification, 
can be applied to any field theory with a broken symmetry, at temperatures below symmetry 
restoration. 

That the dispersion relation of a mode can be expressed fully in terms of static correlation 
functions is nontrivial, but by no means unprecedented. We recall the sound waves, whose 
velocity is u = (dp/ <9e) 1//2 , where p and e are the pressure and the energy density, respectively. 
The sound speed, while being a real-time quantity (the pole in the correlator of the energy 
density T 00 ), can be determined solely from thermodynamics. The relation between the 
speed of sound and the thermodynamic functions is an exact consequence of the existence 
of the hydrodynamic description. A less familiar example is the variational Feynman-Bijl 
formula which relates the phonon spectrum in superfluid helium to the static density- density 
correlation function [[§ . The example most closely related to our problem, however, is that of 
spin waves in antiferromagnets H : the velocity of spin waves at any temperature below phase 
transition is equal to the ratio of the stiffness and the magnetic susceptibility in a direction 
perpendicular to magnetization. Both quantities can be defined from the static response 



Some progress, however, may have been made recently 
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of the system to external fields. The only difference between QCD and antiferromagnets is 
that in the former case the symmetry is SU(2)y x SU(2)^ ~ 0(4), which is spontaneously 
broken to SU(2)y — 0(3), while in ferromagnets 0(3) is broken down to 0(2) || |6|. 

The paper is constructed as follows. In Sec. || we summarize the findings of the paper. 
In Sec. |II| we use simple, but nonrigorous, arguments relying on an effective Lagrangian to 
understand these results. In Sec. [TV] the results are derived in a more rigorous way from a set 
of assumptions about the real-time correlation functions, which comes from hydrodynamics. 
In Sec. [V] we show that our result holds for the simplest field-theoretical model of a scalar 
field theory with broken symmetry. In the Appendix we give a simple derivation of the 
known result about the dynamical critical exponent z, and derive the critical scaling of a 
diffusion coefficient. 



II. SUMMARY OF RESULTS 

We claim that, in QCD with two light flavors, at temperatures below the chiral phase 
transition, the real part of the dispersion relation of sufficiently soft pions is given by the 
following equation: 

u 2 p = u 2 (jp 2 + m 2 ). (2.1) 

In this paper the following terminology is used: u is the pion velocity (although, strictly 
speaking, it is the pion velocity only when m — 0), and m is the pion screening mass (we 
shall show that it is the same screening mass as defined on the lattice). The energy of a 
pion at p = 0, m p = um is called the pion pole mass. 

At finite temperature, the meaning of soft pions may need some clarification. We shall 
understand Eq. fl2.1|) as a statement that the correlators of operators carrying pion quantum 



numbers have a pole at the frequency with the real part determined by Eq. ( |2.1| ). 

The pion velocity u is the ratio of two statically measurable quantities, the temperature- 
dependent pion decay constant f 2 and the axial isospin susceptibility X I5 , 



f 2 

u 2 = ^-. (2.2) 

X I5 



The axial isospin susceptibility X I5 can be defined as the second derivative of the pressure 



with respect to the axial isospin chemical potential (see Sec. |T|). Equivalently, it can be 



defined via the static Euclidean correlator of the axial isospin charge densities, 

l/T 



where 



5 ab X I5 = jdr Jdx (A a (r, x)A b o (0, 0)> , (2.3) 



^^7 7 5 y^, (2.4) 
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ip is the quark field, a,b = 1,2,3, r a are isospin Pauli matrices, Tr T a r b = 25 ab , and (• • •) 
denotes thermal averaging, which can be taken by evaluating a Euclidean Feynman path 
integral with appropriate boundary conditions. The quantity X Ib has been considered pre- 
viously @. 

The pion decay constant / and the screening mass m can be determined from the static 
Euclidean pion correlators at small momenta, which are predicted to have the form 

l/T 

Idr [ dx e-*r"V(T, xy (0, 0)) = ±- 2 ^— 2 , (2.5) 
J J f 2 q z + m z 

o 

where the scalar field ip a is defined as 

f = (2.6) 

Equation (|2.5|) is supposed to be valid when \q\ <C m a , where m a is the (temperature- 
dependent) screening mass of the a meson. 

The way m enters Eq. ( |2.5|) explains why we term it the pion screening mass. As far as 
/ is concerned, the conventional definition of the pion decay constant, f^p^ = (0\A^\n(p)) , 
cannot be used at finite temperature, since neither the vacuum state |0) nor the one-pion 
state \ir) allows generalization in thermal media. Both y/X I5 and / approach f n as T — > 0, so 
both can be viewed as the generalization of f n to finite temperature. We, however, reserve 
the name "pion decay constant" for the quantity / defined in Eq. ( |2.5| ). 2 In contrast to the 
susceptibility X I5 , f cannot be defined at temperatures above critical. 

We shall also show that the Gell-Mann-Oakes-Renner (GOR) relation can be generalized 
to finite temperature, and does in fact become two separate relations for the screening and 
pole masses of the pions, 

f 2 m 2 = X I5 m 2 p = -m 9 (#) . (2.7) 

For simplicity, in this paper we assume m u = rrid = m q . 

Equation (|2.2|) is the direct analogue of a similar equation for the velocity of spin waves 
in a quantum antiferromagnet at temperatures below the phase transition [f|. The coun- 
terpart of / is the stiffness (denoted as p s in Ref. [|J), while X I5 is similar to the magnetic 
susceptibility in a direction perpendicular to magnetization. 

An important point we wish to emphasize is that the relations ( ^.1| )-( P77|) are exact in 
the chiral limit, i.e., the limit when m n (vacuum pion mass) and p are infinitesimally small, 
at any temperature in the interval from to T c . Although the methods we use to derive 
them may look similar to chiral perturbation theory, unlike the chiral perturbation theory, 
we do not consider perturbations around = 0, T = 0, treating both m n and T as small 
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parameters. Rather, we consider perturbations around = at a fixed T, the latter not 
assumed to be small. Therefore, unlike in chiral perturbation theory, we are not able to 
calculate the temperature dependence of parameters such as %■> f ■> u -. or '4 ) '4 ) - However, we 
are able to show that their temperature dependence must be such that relations (|2.1|) - (|2.7|) 
hold. 

Given these equations, exact in the chiral limit, it is legitimate to ask how useful these re- 
lations can be at the physical value of the pion mass m n = 140 MeV. To answer this question 
quantitatively, a (perhaps lattice) calculation of the T and dependence is needed, which 
is beyond the scope of this paper. However, we shall attempt to give a semiquantitative 
answer to this question. 

It is easier to begin with T = 0. In this case formula (|2.1| ) becomes trivial: the dispersion 
relation is exact, simply by Lorentz invariance, for any p and m = m n . It is a straightforward 
exercise to show, using PCAC (partial conservation of axial vector current), that to leading 
order, i.e., C(m^), the correlator in Eq. ( [2. 3D whose value at m n = is X I5 , does not depend 
on m n . Thus one can expect that measuring the correlator ( p. 3D even at the physical m n 
(or higher, as is typical in a lattice calculation due to the high price of simulating with 
light quark masses) is not significantly different from X I5 . In order to extract the value 
of / from the relation (|2.5| ) reliably, one requires that higher mass states in the channel 
with pion quantum numbers do not significantly contaminate the exponential falloff of the 
correlator in coordinate space. Since the masses of such states are typically above 1 GeV, 
the exponential tail contribution from the light state of mass m n = 140 MeV should be easy 
to separate. 

At nonzero T the region of applicability of Eq. (|2.2j ) is limited by terms of higher powers 
in p. For simplicity, in the exact chiral limit m n = the dispersion relation takes the form 
u = \p\ — jD'p 2 + • • • , where D' is a temperature-dependent parameter (see Sec. [TV] and 
the Appendix). The condition on momenta which is required to neglect nonlinearity in the 
dispersion relation is \p\ <C 1/D'. Calculation of the diffusion constant D' is a challenging 
task. Taking an estimate from Ref. ||, D' — CT 3 /f£, with a numerically rather small 
C ~ 0.1, we can conclude that even for T ~ f n at momenta of order \p\ ~ 100 MeV the 
nonlinearity in the dispersion relation is still small. 

As T = T c the dispersion relation is essentially nonlinear, which manifests itself in the 
divergence of D' as T — > T c (see the Appendix). As T — > T c the maximum momentum 
at which the dispersion relation can be considered linear (in the chiral limit) decreases 
and vanishes at T c . Although the power with which the width of the linearity window 
shrinks to zero can be determined (by extension of arguments given in the Appendix), the 
pre-exponent is unknown, and a nonperturbative calculation is, in principle, required to 
determine quantitatively the size of the nonlinearity in the dispersion relation at a given T 
and |p|. Such a calculation is beyond the scope of the paper. 

Equations (|2.1|) - ([2.7|) were used in Ref. || to extract information about the critical be- 
havior of the pion velocity and masses near the critical temperature. A brief derivation of 
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Eqs. ( |2.1|) - (|2.7|) was also sketched in Ref. 0. We present a more extended version of this 
derivation in Sec. DTD. 

In Sec. [IV| we provide a new derivation of the relations (|2.1|) -( [2.7|) using the operator 
approach based on hydrodynamic equations. For this purpose we shall need expressions 
for X I5 and / in terms of equal-time rather than static (zero frequency) correlators. These 
relations follow from definitions ( |2.3| ) and fl2.5p in the exact chiral limit (m n = 0) or when 
the temperature is high enough (T 3> m n ): 

5 ab X I5 = ^ jdx(A a (t,x)A b (t,0)). (2.8) 

In the chiral limit, the equality of the static correlator [Eq. ( |2.3| )1 and the equal-time cor- 
relator [Eq. ( |2.8| )1 of Aq is a consequence of the conservation of the axial isospin charge. 
Similarly, 

I fdx e-**> B (f , x)ip% 0)> = ±- 2 -^— 2 . (2.9) 
1 J j Q + m 

Note that, in contrast to Eqs. ( |2.3| ) and (pT5|) , Eqs. ( |2.8| ) and ( p.9| ) require an additional 
condition T 3> m n . This condition is, however, not so dramatic, if one recalls that the 
relavant energy scale of order T is the lowest nonzero Matsubara frequency, InT . Thus, in 
practice, one requires 27rT ^> m 7T , which is satisfied reasonably well for T of order T c « 160 
MeV and the physical pion mass. 

However, a much stronger condition is required for the applicability of the hydrodynamic 
description: \p\ (and m n ) must be smaller than the typical collision rate r _1 ~ T 5 /f£ 
(according to Ref. WW). At small T <C fn this condition is much stronger than the conditions 
that are required for the effective Lagrangian derivation to hold (typically \p\,m n <C m a - 
see Sec. |TTT|) . However, both derivations apply in the required regime infinitesimally close 
to the chiral limit. The fact that the hydrodynamic derivation has a smaller validity range, 
however, does not mean that the result is not valid outside this range. As an example, 
consider the case of T = 0, where the results (|2.1| ), (|2.2|) hold trivially, while hydrodynamics 
does not apply at all. A less trivial example is presented in Sec. 0, where explicit calculation 
verifies Eq. ( |2.2|) . This weak-coupling calculation, however, does not require any notion of 
hydrodynamics. 

The advantage of the hydrodynamic approach is that it allows us to consider properly 
the effects of the dissipation, which are inherently beyond the Lagrangian approach. Before 
proceeding to the new operator derivation (Sec. |TV| ) we review the derivation based on the 
effective Lagrangian approach. 
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III. EFFECTIVE LAGRANGIAN APPROACH 



A. Axial isospin susceptibility and f n at T = 

Since an important role in our analysis is played by the axial isospin susceptibility, we 
first consider this quantity at zero temperature and its relation to the pion decay constant. 
The quark part of the QCD Lagrangian at finite axial isospin chemical potential fi I5 is given 
by 

_ 7g 

£quark = i^D^ ~ m q 1p1p + £4/5^7075 7^0 , ( 3 -l) 

where is the color covariant derivative. This chemical potential /i /5 is coupled to the axial 
isospin charge [cf. Eq. fl2~lD l, which generates the SXJ(2)a part of the SU(2)y x SU(2)a 
chiral symmetry. 

The response of the QCD vacuum to fi I5 can be found from the effective chiral Lagrangian. 
The latter, to lowest order of momenta, masses, and chemical potential, is completely fixed 
by the chiral symmetries and is given by 

C cS = ^TrV v SV v S t + ^-^ReTrS, (3.2) 

where E is an SU(2) matrix whose phases describe the pions, E = e lTa7r °'^ 7r , and V denotes 
the covariant derivative, which is defined as 

V E = «9 E - ^ /5 (r 3 E + Er 3 ) , V*E = <%E , t = 1,2,3. (3.3) 



At /j, I5 = 0, the Lagrangian ( |3.2| ) is the standard chiral Lagrangian with two phenomeno- 
logically determined constants f n and m^. The way /x /5 enters the effective description is 
completely fixed by symmetries to lowest order. This can be seen by promoting the SU(2)a 
symmetry to a local symmetry and treating /j, I5 as the time component of the SU(2)^ vector 
potential | fH| . 

The Lagrangian (|3.2| ) and its derivation is analogous to the case of the effective Lagrangian 
at finite (vector) isospin chemical potential fij considered in Ref. |TT[] . A significant difference 



between the cases studied here and in Ref. |Tl| is that the QCD vacuum breaks the SU(2)^ 
(axial isospin) symmetry spontaneously, but remains symmetric under the SU(2)y (vector 
isospin) symmetry. It is important to note, however, that the SU(2)a is a symmetry of the 
Lagrangian (at m q = 0), as good as the SU(2)y. In particular, the axial isospin current A® 
is conserved in the chiral limit. Thus, it is entirely legitimate to consider the theory at finite 
fM I5 and use symmetry arguments to fix the fM I5 dependence of the effective Lagrangian. 

The vacuum energy density depends nontrivially on /x J5 already for arbitrarily small /i J5 . 
(This is in contrast to the case of the isospin chemical potential where [i I needs to be 
larger than a threshold equal to m n in order to change the ground state.) The isospin axial 
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susceptibility, at fi I5 = 0, is easy to determine using the effective Lagrangian Q3.2|) : 



— 



d 2 £ 



15 



M/5 = 



d 2 C 



off 



J 7T 



(3.4) 



I?, 



S=l 



This result is of potential importance for lattice QCD calculations because, in principle, 
it allows one to determine f n directly by measuring the axial isospin susceptibility. To our 
knowledge, this has not been done on the lattice for temperatures below the chiral phase 
transition. The isospin susceptibility (as well as the isoscalar, i.e., baryon number suscep- 
tibility) measurement has been done using staggered fermions in |[L2| . In the formulation 
used in Ref. [[12|], introducing the axial isospin chemical potential would require replacing 
exp(— jia) on the links (where a is the lattice spacing) with exp(— fia( x ), where ( x is the usual 
staggered factor ( x = (— l) x 1+X2 +X3+X4 — the representation of 75 in the staggered fermion 
action. 



B. Pion velocity 



To obtain Eq. ( [2.2| ), we expand the previous discussion to nonzero temperature. We 
first presume that the dynamics of the pions is described by some effective Lagrangian 
£ e ff. Strictly speaking, this is not correct since dissipative effects cannot be included in 
the effective Lagrangian. We can expect to recover the correct answers if, in the infrared, 
the pion thermal width is negligible compared to its energy. This has been seen in explicit 
calculations at low T, where chiral perturbation theory can be used [13, |14|]. We furthermore 



assume that this Lagrangian is local and can be expanded in powers of momenta. To lowest 
order, the Lagrangian is fixed by symmetries up to three coefficients, ft, f s , and m, 

C cS = IJ-T T V S VoS 1 " - ^Tr d l Y>d l ^ + — ^Re Tr S . (3.5) 

Due to the lack of Lorentz invariance, f 2 and f 2 are independent parameters. The covariant 
derivative Vo is the same as defined in Eq. (|3.3|) . The dispersion relation following from this 
Lagrangian has the form ( p.l| ) where 

u = fy. (3.6) 
ft 

At zero temperature ft = f s = fir, and the pion velocity u is equal to the speed of light. We 
now show that, at finite temperature, all three parameters ft-, f s , and m can be determined 
from equal-time (or static) correlation functions. 



Repeating the same argument as in Sec. |III A| , we can show that ft is related to the axial 
isospin susceptibility, 

X I5 = fl , (3.7) 
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but now X I5 is denned as the susceptibility at finite temperature. Note that, as a suscepti- 
bility with respect to a conserved charge, X I5 is free of ultraviolet divergences. 

We now need a prescription to compute f s and m from static correlation functions. In 
order to make the connection, we generalize the mass term in Eq. (|3.1|), 



£quark = ty-fD^ ~ {lp L M^ R + H.C.) + fi^A* , (3.8) 

and regard M as an external field, M = M{x). The Lagrangian (|3.8|) possesses the symmetry 
ifti, -> LipL, ?p R -> RipR-i M -> LMR\ where L, R E SU(2). If M(x) is a slowly varying 
function of x, its effect can be captured in the effective chiral Lagrangian. The requirement 
that the effective Lagrangian preserves this symmetry fixes the form of its mass term, 

P f 2 1 — 

£ cfr = ^-Tr VoSVoS 1 " - ^Tr d^d^ - - (^)Re Tr M f S . (3.9) 

We shall limit ourselves to a particular ansatz of the external field M(x): 

M(x) =m q e taa{x)ra . (3.10) 

The second derivative of the partition function with respect to a a can be computed in 
both the microscopic theory ( |3.8| ) and in the effective theory ( p.9| ). In the microscopic theory, 
we find 

where ip a is defined in Eq. (|2.6|) . On the other hand, from the the effective Lagrangian ( |3.2| ) 
we find 

5a«(x)5a b (0) = m ^ )2 ' where f( x ) = ReTr ^(x)/2 . (3.12) 

Comparing Eqs. ( |3.11| ) and ( |3.12j ), we see that correlation functions of ip a (x) defined in the 
microscopic theory and of <p a (x) defined in the effective theory are equal. This equality 
should hold for small momenta when the effective theory is applicable. 

The correlation function of a (x), on the other hand, can be calculated by expanding the 
effective Lagrangian to second order in <p a . The Matsubara propagator of <p a is 

l/T 

/Sjab 
dxe^{r{x)m) = f 2 q 2 + mq2 + m2 y 9b = ^Tn. (3.13) 
o s 

This propagator, for q = and small q, should be equal to the propagator of (p a . This 
establishes Eq. Q2.5Q , with the identification / = f s . Together with Eqs (|3.6| ) and (p.7|), 
this is our result for the pion dispersion relation. Furthermore, it is natural to assume that, 
for small q, the dynamics of <p a is slow, so at high enough temperature (T ^> m n ) one can 
regard (p a (r, x) as independent of r. In this case the left hand side of Eq. (|3.13|) is equal to 
that of Eq. (pi). 
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It is instructive to write Eq. ( |2.9| ) in coordinate space, 



(y B (0,gV(0,O))= 6 {\x\^T-\m- x ). (3.14) 

We see that measuring, in the microscopic theory, the large-distance equal-time correlation 
function of the tp a defined in Eq. Q2.6Q we can extract the screening mass m and the decay 
constant /. Combined with the determination of the susceptibility X I5 and Eq. ( |3.7|) , the 
dispersion relation of soft pions is now completely known. 

Below we shall provide a more systematic proof of these relationships using an operator 
approach, with crucial inputs from the hydrodynamic theory. We shall also demonstrate the 



validity of the relation ( |2.2|) in an explicit lowest order perturbative calculation in the linear 
sigma model. 



IV. HYDRODYNAMIC (OPERATOR) APPROACH 

From a modern perspective, hydrodynamics is an effective theory operating at sufficiently 
large distance and time scales. (By "sufficiently large" normally we mean scales larger than 
the mean free path, or the relaxation time.) As such, it is the most suitable framework to 
discuss low-energy degrees of freedom (like pions) at finite temperatures. This modern point 
of view, as opposed to the view of hydrodynamics as a purely phenomenological description, 
has been in existence for a long time J0| . From this philosophy, it is not surprising that the 
same hydrodynamic theory describes systems with very different microscopic dynamics, and 
that systems with different symmetries (like normal fluids and superfluids) correspond to 
different hydrodynamic theories. Similarly to the effective Lagrangian, hydrodynamic equa- 
tions can also be viewed as a particular way to satisfy Ward identities at finite temperatures 
T6fl . In our case, hydrodynamics is nontrivial due to the chiral symmetry breaking [17]. We 



will derive the constraints placed by hydrodynamics on the dynamics in our problem, and 
show how our results on the pion dispersion relation follow from there. Our treatment is 
similar, to some degree, to that of Ref. [|J. 



A. Basic assumptions of hydrodynamics 

We shall assume that, as one goes sufficiently far into the infrared, the dynamics of 
any interacting finite-temperature system can be described in terms of a finite number of 
fields, which will be called hydrodynamic variables. To be relevant in the infrared, the 
fluctuations of these fields (either thermal fluctuations, or those due to external sources) 
should relax arbitrarily slowly. This requirement eliminates most of the degrees of freedom, 
which typically relax during a relaxation time determined by the microscopic dynamics. 
However, the following fields are obvious hydrodynamic variables: 
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(i) The densities of conserved quantities, including the energy density T 00 , the momentum 
density T° J , and the densities of conserved global charges (i.e., the zeroth components 
of conserved currents). These fields cannot relax quickly because of the conservation 
laws. A configuration where charges fluctuate over a length scale L much larger than 
the mean free path can relax only by diffusion, which takes place over a time propor- 
tional to L 2 . The relaxation time diverges with the wavelength of the perturbation. 

(ii) The phases of the condensates which break global symmetries. At zero temperature 
the fluctuations of such phases correspond to Goldstone bosons, whose energy can be 
arbitrarily small. At finite temperature below symmetry restoration, one should also 
expect the long-wavelength fluctuations of the condensate phases to relax slowly. 

(iii) Near the critical temperatures of second-order phase transitions, the order parameters 
themselves (not just the phases) should be considered hydrodynamic variables 



For example, it is believed that, in QCD with two massless flavors, the chiral phase 
transition is of the second order, where the a meson becomes degenerate with the pions. 
Near T c , hence, a should be included in the hydrodynamic description. In contrast to 
the fields in the categories (i) and (ii), the rate of relaxation of order parameters is 
controlled by the closeness of T to T c , but not by the wavelength of perturbations. 3 

If one makes a further assumption that the fields listed in (i)-(iii) exhaust all slowly 
relaxing ones, then the set of hydrodynamic degrees of freedom is completely fixed once all 
symmetries of the theory and the pattern of symmetry breaking at the given temperature 
are known. In this paper we shall limit ourselves to temperatures far away from any second 
order phase transition, so order parameters are excluded from hydrodynamics. This set of 
hydrodynamic variables always contains the energy and momentum densities T 00 and T°\ 
For QCD below the chiral phase transition, one has, in addition, the densities of baryon, 
isospin, and axial isospin charges, and the phases of the chiral condensate. 

After identifying the hydrodynamic degrees of freedom, one can proceed in various ways. 
One can ask about the the equations of motion that the hydrodynamic variables obey. For 
our case, it is a rather nontrivial task, because of the multitude of the fields involved. In 



Ref. ||17|| , the Poisson-bracket technique is used to derive the dissipationless hydrodynamic 
equations. One finds a system of fully nonlinear coupled differential equations, generalizing 
the equations of fluid dynamics and equations of motion of the nonlinear sigma model. 
However, the connection of this procedure to the fundamental (microscopic) field theory 
has not been made, and the physical meaning of several temperature-dependent parameters 



appearing in the final equations (called f t , f s , and f v in Ref. lfL7| ) is not at all clear 



3 We shall not consider the possibility of Abelian gauge fields in the Coulomb phase (as opposed to the 
Higgs phase) present in magnetohydrodynamics. 
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Alternatively, one can ask the question: what are the constraints that hydrodynamics 
places on the correlation functions of the hydrodynamic variables? It is clear from the discus- 
sion above that the real-time correlators of the hydrodynamic variables [i.e., (0(t, x)O(0, 0)), 
where O is the variable] are long range, i.e., have power-law (but not exponential) falloff, 
at least in the timelike regime t —>■ oo, x = fixed. 4 It is less trivial to decide about the 
correlators of the operators not belonging to the set of hydrodynamic variables. To this end, 
the following additional assumption is made. 

(iv) All local operators can be expressed as local functions of the hydrodynamic operators 
and their spatial derivatives, up to corrections which have short-ranged correlations 
that go to zero exponentially when either temporal or spatial separation goes to in- 
finity (t or x —>■ oo). (These short-range parts correspond to the "noises" in the 
hydrodynamic theory.) 

Physically, these assumptions mean that the dynamics in the infrared can be described in 
terms of hydrodynamic variables only, and is equivalent to the assumption of local ther- 
modynamic equilibrium: the values of all variables are determined by specifying a few. 
The assumptions (i)-(iv) form the starting point of our construction of the hydrodynamic 
equations. 

B. Linearized hydrodynamics for soft pions 

To be less abstract, let us consider QCD in the chiral limit, at a finite temperature below 
the chiral phase transitions, and zero chemical potentials. The hydrodynamic operators are 
the energy density T 00 , the momentum density T°% the baryon density ip^ip, the densities 
of vector and axial isospin charges, 

^ a = ^7°y^, ^ = ?7°7 5 y^, (4.1) 
and the pion field, defined as Eq. (|2~6|) , 

<P a = • (4-2) 

As we shall see below, to linear order the dynamics of Aq and (p a is decoupled from other 
modes. 

Let us consider the equation for Aq. Infinitesimally close to the chiral limit, we can derive 
from the QCD Lagrangian the familiar PCAC relation: 

d^ = m q (^W a - (4-3) 

4 The hydrodynamic correlators are not required to have power-law decay in the regime t = fixed, x — > oo. 
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The left hand side contains, in addition to Aq, the spatial components of the axial current 
A ai , which are not hydrodynamic variables. According to the assumption (iv), we can 
express A at as local functions of the hydrodynamic operators and their spatial derivatives, 
plus a short-ranged part. If we work to leading order in the power of the fields (which means 
the linear order), the only ones suitable are (p a and Aq which are parity odd and isovectors. 
The spatial index in A ai forces one to have at least one spatial derivative. Thus, to leading 
orders in fields and derivatives, the only terms consistent with symmetries are 

A ai = _fQ i(p a _ Dd , A a _ ^ai ? ^ ^ 

where f 2 and D are coefficients depending on the temperature, and £f is the short-range 
("noise") part of A m . (At this step, we have not yet related f 2 to the static correlation 
functions. We will do it later on.) Equation (|4.3|) now takes the form 

d A a = fV V + mS^ a + DV 2 A a , + d t C ■ (4.5) 

The parameter D can be interpreted as the diffusion coefficient for Aq] however, Eq. ( |4.7a| ) 
is more complex than a diffusion equation. Now let us discuss the equation for ip a . The time 
derivative of (p a , not being a hydrodynamic operator, can be expanded as 

d oV a = —Aq - n' lV a + k 2 V V + V a , (4-6) 

where X, k[, and k 2 are again coefficients dependent on temperature, and rj is a short-range 
noise. We have kept the term with the second derivative of ip a as well as the term with no 
derivative. The reason for doing so is that k' x is suppressed by the quark masses. Indeed, 
in the chiral limit the state with ((p a ) ^ 0, (Aq) = can be another vacuum which stays 
unchanged with time. This is consistent with Eq. ( |4.7b| ) only if k' x = in this limit, thus at 
small quark masses k[ is small. For this reason we keep the V 2 ip a term. 

Introducing the parameter m 2 defined so that f 2 m 2 = —m q (%p%l)) (which we still have to 
relate to static correlators), and Ki so that k[ = Kim 2 , Eqs. ( |4 . 5| ) and ( |4.6| ) can be written 
as 

d A a = / 2 (V 2 - m 2 )ip a + D\7 2 A a Q + d£ ai , (4.7a) 
d p a = -A a Q + (k 2 V 2 - K im 2 ) V a + rf . (4.7b) 

Equations ( |4.7| ) are the linearized hydrodynamic equations governing the evolution of Aq 
and ip a . The correlation functions of Aq and ip a can be found if one knows the correlators of 
the "noises" £ a * and rf. By construction, these fields have only short-range correlations, so, 
if one is interested only in the dynamics at large distance and time scales, these correlations 
can be replaced by delta functions. Isospin symmetry and rotational invariance require the 
correlators to be of the following forms: 

(C(-r)i h! i!))) = F^5\x - y) , (4.8a) 
(V a (x)v b (y)) = F v 5 ab 5\x - y) , (4.8b) 
(r(x)r, b (y)) = 0. (4.8c) 
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up to corrections proportional to derivatives of 5 4 (x — y) which will be neglected since they 
are of higher order in momentum. Equations ( J4.7| ) and fl4.8[ ) completely determine the 
hydrodynamics of soft pions, to the linearized order. 



C. Relation to static correlators 



Our next task is to relate, as much as possible, the parameters appearing in Eqs. 
and ( |4.8| ) with the equal-time correlators of Aq and ip a . 

First, we multiply Eq. ( |4.7b| ) by A b , taken at the same time moment, and integrate over 
space. One finds 



/ 



dx(ip a (t,x)A b ( y t, 0)) = jdx 



l(A a (t, x)A b (t, 0)> - K im 2 ( V a (t, x)A b (t, 0)> 



r uw ' uw " u (4.9) 

+( v a (t,x)A b (t,o))] . 



In the right hand side, the second term is proportional to m 2 , which is small in the chiral 
limit, and hence can be neglected. The last term in the integrand will be shown later to 
vanish, 

(r l a (t,x)A b (t,0))=0. (4.10) 

Equation ( |4.10| ) is trivial if understood as (r] a (t + e, x)A b {t ) 0)) = 0: A b Q cannot be correlated 
with the noise r\ in the future. What is somewhat less trivial (and will be checked a posteriori) 
is that (r) a (t — e, x)A b ) (t, 0)) also vanishes: the equal-time correlator of r\ and A does not 
depend on how the equal-time limit is taken. Therefore, only the first term survives; and 
by definition of X I5 it is equal to 

Jdx^(A a (t,x)A b (t,0)) = ^Y 1 - (4-11) 

On the other hand, the left hand side of Eq. (4^) can be computed explicitly. To this 
end, we write 

(^x)A b (0))=iTT{e^ H [H,^(x)}A b (0)}. (4.12) 

(we drop the time variable t which is an argument of all operators). We now make use of 
the following expansion: 

[e-^, V a ] = -Pe-f H [H, V a ] - ^e-? H [H, [H, ^]] + ■ ■ • . (4.13) 

The expansion parameter here is f3qo = qo/T where go is the frequency of variation of (p a . 
Since we are dealing with the low-frequency modes in (p a , we can ignore all terms beyond 
the first in Eq. ( f4.13|) . Equation ( [4.1 2|) now becomes 

m*)<(0)) = -iTTr{[e-P H ,<p a (x)}A b (0)} = ^([^(x), A b (0)]) = T5 ab 5*(x) . (4.14) 
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In the last transformation, we make use of the commutation relation 

[^(x),A b (0)}=z6 ab 6 3 (x)p-. (4.15) 
Comparing to Eq. ( |4.11|) , we find 

X = X I5 . (4.16) 

Thus, we show that X is the axial isospin susceptibility of the system, which is a static 
quantity. The proof we have just presented is similar to that of the equipartition theorem 
in statistical mechanics. 

Analogously, we can show that 

J dx e- iq - x (A a (x)cp\0)) = -T5 ab . (4.17) 

On the other hand, from Eq. (4.7a) 



I dx e- iq - x (A*(x)ip b (0)) = - f{q 2 + m 2 ) ( dx e- iqx (ip a (x)ip b (0)) 
J J (4.18) 

- Dq 2 / dxe~ iq x (A a Q (x)ip b (0)) , 

where we dropped the term proportional to 

aV (*,())) (4.19) 

as it will be shown to vanish in the same manner as (r] a A b ) inEq. (|J). Moreover, Dq 2 (Afa b ) 
can be neglected compared to {Aq^) ~ q (AQ(p b ), since we expect the pion to have a linear 
dispersion relation and, for small enough momenta, q ^> Dq 2 . Equating Eqs. ( |4.17| ) and 
Q4.18] ), we find that equal-time correlators of the pion field must have the form of a Yukawa 
potential, 

/T l) ab 
dx e-^^y (0)> = j 2 ^^ 2 ■ (4-20) 

Thus, the parameters j 2 and m 2 appearing in Eq. (f4.7a|) are the same ones defined in Eq. 
( |2.9| ) via the equal-time correlator of (p a . 



The parameters D, K\, k 2 , and F n cannot be expressed individually in terms of the 
equal-time correlators, but some relations between them will be derived below. 

D. Hydrodynamic correlation functions 

From Eqs. ( |4.7|) and ( [4.8|) one can easily compute the real-time correlators of A$ and (p a : 

#xe** (A»(x)A b (0)) = ^- q M^W^\ I iF2l > ( 4 ' 21a ) 

[(g - uJ q ) 2 + ir 2 ][(g + uJq) 2 + l^q] 

d'x (<p a (x)<p\0)) = 5 ab - , , (4.21b) 

l(qo - u q ) 2 + z^ll [(go + uj q ) 2 + i r2] 
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where 



u\ = ^(q 2 + m 2 ), (4.22a) 
T q = K im 2 + (D + K 2 )q 2 . (4.22b) 



The correlation functions (4.21) peak around go ~ i^gj an d the width of the peaks r q is 
much smaller than uo q at small q. The correlators have poles corresponding to the pion 
collective excitations with the dispersion relation q = uj q — ^T q . We note, moreover, that 
all parameters entering the real part of the dispersion relation ( |4.22a| ) can be determined 
from the static correlation functions. 

It is now possible to explicitly check that the noise-field correlators appearing in Eqs. 
Q4.10Q and (|4.19|) vanish. Indeed, by using Eqs. tfL7§ and (|4~8| ) we find 

<„(*, M 0)) = **/ (gj e-<-' + - , (4 . 23a) 

<«.(«, * M t>, 0)) = -d-/ ^ >+- • (4.23b) 

The integrands in Eqs. fl4.23|) have two poles, both located in the upper half plane: go = 
±uj q + ^T q . When t > t', when taking integrals over g , one can close the contour in the lower 
half plane, so the integrals are obviously zero, which corresponds to our previous remark 
that the fields cannot correlate with the noises in the future. If one takes t —>■ t' from below, 
the integrals are also zero since they are equal to the sums of the residues of the integrands 
(which are zero because these functions behave as q^ 2 at large go.) Therefore, we have 
checked that the equal-time noise-field correlators ( 4.10|) and ( f4.19| ) indeed vanish. 



One can find the relations between the amplitude of the noise correlators, F% and F v , and 
the parameters characterizing the damping in Eq. fl4.22b| ). One integrates Eqs. ( |4.21|) over 
go and obtains the equal-time correlation functions 

,2 77 i v 2 , ,2. 



J dx e-^ x (A a (t, x)A b (t, 0)> = 5 ab H I ^ q v , (4.24a) 

/ dx e-^{<p a (t, x)<p% 0)> = S ab q2F ^ 2 X \ 5 f qFv . (4.24b) 

Comparing these correlators with Eqs. ( |2.8| ) and (|2T9|), taking into account Eqs. ( |4.22|) , one 
finds 



Ft = 2TX I5 (D + k 2 - «i) , (4.25a) 
2Tk 1 

1 

Our results for the hydrodynamic correlations functions are given by Eqs. (|4.21|) and ( [4.251) . 



For completeness, we give here also the result for the cross correlator of ip a and Aq\ 

d'x e<** ( V a (x)A b (0)) = , 2 J^f\ ,2^i F 2i • ( 426 ) 

[(go - ujq) 2 + [(g + u q ) 2 + z r q \ 
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It is instructive to compare the correlator of the axial isospin charge density Aq [Eq. 
( |4.21a| )] with that of the vector isospin charge density V a . The dynamics of the latter is 
completely diffusive and is given by the equations 

d V Q a - DjV 2 V a = diC ai , (4.27a) 
(C(x)C bj (0)) = 2TD I X I 5 ab 5 ij 5\x) , (4.27b) 

where D 1 is the diffusion constant for the isospin charge, and Xj is the (vector) isospin 
susceptibility. The correlator of V a is 

/9TD y n 2 
d*xe^(V o a (x)V o \0)) = . (4.28) 

% u iH 

The pole of the correlator is located at purely imaginary frequencies go = ^i-Djq 2 , as it 
should be for a purely diffusive mode. 

Using Eqs. ( |4.28|) and Q4.24|) we can understand the limitations of the Lagrangian ap- 
proach and the role of the dissipative processes. For example, one could try to apply the 
method of Sec. HJ to determine the vector isospin susceptibility Xj, introducing an isospin 



chemical potential /i/. At zero temperature such a method was used in Ref. ||11|| . As ex- 
pected, the isospin susceptibility Xj vanishes at T = since it takes finite energy to excite 
isospin degrees of freedom (pions) and change the isospin density. However, at nonzero 
temperature, the naive application of the effective Lagrangian method would predict that 
Xj — also. 5 On the other hand, we should expect Xj ^ at finite temperature, even 
at very small T, because of the presence of isospin-carrying pions in the thermal medium. 
Looking at the correlator ( 4. 28] ), one sees that the equal-time correlator defining Xj, being 



the integral of Eq. ( |4.28|) over go m the limit q — > 0, receives the main contribution from 
very slow (diffusive) modes: go ~ Djq 2 . However, the effective Lagrangian (|3.5| ) only de- 
scribes faster (propagating) modes of the if field with g = uj q ^> D t q 2 . The slow diffusive 
modes which contribute to Xj are not present in the Lagrangian ( |3.5| ). For the axial isospin 
susceptibility X I5 , the situation is different, for the integral of Eq. Q4.28D over frequencies in 



the limit q — » is concentrated entirely near values go = ±lj q . Thus, unlike Xj, the value of 
X I5 at finite T can be found correctly using the effective Lagrangian approach. 

E. Comparison to previous results 

As a by-product of our analysis, we obtain a generalization of the Gell-Mann-Oakes- 
Renner relation to finite temperature: 

f 2 m 2 = -m q l^ip) . (4.29) 



We thank T. Schafer for pointing this out to us. 



17 



This equation has the same form as at zero temperature. At finite temperature, m should be 
understood as the pion screening mass, and the exact meaning of the temperature-dependent 
"pion decay constant" f 2 is given by Eq. (|4.20 ). One can also write the GOR relation in an 
alternative form, 

X I5 m 2 p = -mS^) , (4-30) 

where m p = um is the "pole mass" of pions. 

That the velocity of pions at finite temperature is different from the velocity of light has 
been seen in the second order of chiral perturbation theory (i.e., in order T 4 / f^) ]TJ| ^D[. The 
authors of Ref. fl]| also introduced two pion decay constants /* and f s , which correspond 
to our ^X I5 and /, and checked the validity of the GOR relation ( |4.30| ). However, these 
constants were defined only in chiral perturbation theory, and only at small T; there has 
not been any attempt to give a precise definition of the constants /* and f s at temperatures 
comparable to T c . We have, in contrast, given a precise meaning to the constants X I5 and 
/ in terms of equal-time correlation functions that can be measured on the lattice. The 
GOR relation now contains only well-defined quantities. Although Eq. ( }4.30|) cannot be 
checked on the lattice, since it contains the pion pole mass, the version (|4.29|) can be verified 
numerically since all quantities entering it are statically measurable. 

We also note that Eqs. (|4.7|) , without the noise and dissipation terms, can be obtained 
by linearizing the hydrodynamic equations obtained in Ref. JI^ by the Poisson-bracket 
technique. One can then identify the parameters f t and f s in Ref. fl7| as f 2 = X I5 , f 2 = f 2 
(the parameter f 2 of Ref. JL7|) is equal to the vector isospin susceptibility). We did not try 
to reproduce the full nonlinear hydrodynamic equation of Ref. |]17fl in the present approach. 



V. AN EXPLICIT EXAMPLE: A LINEAR SIGMA MODEL 



It is instructive to explicitly verify our relation ( |2.2| ) between the velocity of the Goldstone 
bosons, the axial isospin susceptibility, and the temperature-dependent decay constant in a 
model where weak-coupling calculations are possible. The simplest theory with spontaneous 
breaking of a continuous symmetry is the linear sigma model. We shall show that Eq. ( f2.2|) 
indeed holds in this model to leading order of perturbation theory 

It is important to emphasize that the sigma model considered in this section is only 
meant to serve as an example of a theory where we can explicitly check our results by 
perturbative calculations. These models do not describe QCD at T of order T c — a theory 
with no apparent small parameter. We claim, however, that the connection between the 
dispersion relation of the pions and static correlation functions is model independent, and 



can be derived from only a few general assumptions stated in Sec. [TV 
We start with the following Lagrangian 



c = l -d^ a d,r + |W - , (5.1 
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where a = 1,2, ... ,N. We assume A <C 1, so perturbation theory can be applied. When 
/i 2 > 0, the O(N) symmetry is spontaneously broken at zero temperature. We choose the 
vacuum to align in the <p N direction. At zero temperature, and at tree level, the vacuum 
expectation value of <p N is 

(<j> N ) =v = ^-. (5.2) 

At finite temperature the expectation value of <p N is different from t> . Denoting it as v and 
replacing (p N = v + o in the Lagrangian (|5.1|), we obtain 



C = \{d,af + \(d^f - (Xv 2 -^ 2 )va - W-/i> 2 - \(^ 2 -^ 2 

AAA ( 5 - 3 ) 
- Xva 3 - Xvaiz 2 - V - -(tt 2 ) 2 - -a 2 7r 2 , 

4 4 V ; 2 



where 7r = (<p l , <f) 2 , . . . , 7V_1 ). From Eq. (|5.3|) the Feynman rules can be easily written down. 

The value of v is determined by the condition that (a) = 0. For computation, we will use 
the Matsubara (Euclidean) formalism, in which this condition reads, to one-loop order, 

{Xv > _ M> + XvT £ J*_ (_2_ + ^) = o , (5 . 4) 

In this Section, P 2 = p^ + p 2 , p = \p\. The sum-integral in Eq. (5.4) is ultraviolet divergent. 
This temperature-independent divergence can be absorbed into the redefinition of /i 2 . The 
temperature dependence of v comes from the thermal part of the Lagrangian. We shall be 
interested in temperatures of order t> , so T > vn a . In this case, 

Aw 2 - /x 2 + — ^AT 2 = , (5.5) 

or 

T 2 \ ^ u ^ 12 n' 



Equation (|5.6|) is valid everywhere except for a narrow Ginzburg region near T c . 



The masses of a and it are computed, e.g., in Ref. [21|. At the one-loop level m a receives 
contributions from four diagrams 



(where a propagators are denoted by solid lines, and tt propagators are drawn as dashed 
lines). The last two bubble diagrams are negligible in the finite-temperature regime we 
consider, so 



ml = 3A^ - + XvT £ ( J^—t + V 1 I = 2X " 2 ' (57) 
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and thus decreases with temperature. The Goldstone boson receives corrections to its mass 
from three one-loop graphs, 



and remains massless since 



ml = Xv 2 fx 2 - /i 2 + XvT V / ( m 1 - + 



N -1 



Pa 



PO 

dp 



(5.8) 



0. 



mi 



To find the dispersion relation of the pion near q = 0, one has to expand the pion self- 
energy £(g) in powers of go an d <?• The only one- loop diagram that depends on q is the 
bubble diagram, which can be evaluated either by using the Schwinger-Keldysh real-time 
formalism, or by doing calculations in the Matsubara formalism and then performing an 
analytic continuation from imaginary to real go- By either method one finds 



s(g) 



-4AV 



dp 



+[f(co p+q )-m p )} 



[f(uj p+q )+f(Q p )+l} 
1 

+ 



1 



Qp+u p +q— go 



+ 



1 



^p+^p+q+go 



VL p —u p+ q—qQ VL p — cUp+q+go 



(5.9) 



where u p and Q p are the energies of the 7r and a particles with momentum p: uo p = p, 
^p — (p 2 + m l-) 1 ^ 2 j an d /(^) is the Bose-Einstein distribution function f(uj) ■ 



(3u> 



We need to compute the coefficients of q 2 and g 2 in the expansion of E(Q). To compute 
the coefficient of g 2 , one can set go = and Eq. ( p^9|) becomes 



E(0,g) =4AV 



dp 



UJ 



P+Q 



UJ. 



P+Q 



(5.10) 



The right hand side of Eq. ( |5.10| ) can be expanded in powers of q. The constant is compen- 
sated by the other diagrams. The g 2 term gives rise to an integral which is dominated in 
the infrared region p ~ m a , and is of order 



S(0,q) 



X 2 v 2 T 



m°„ 



q 2 = 0(\V 2 )q 2 + 



(5.11) 



In the last equation we assumed T ~ t>o, and m 2 ~ AT, which is valid when T is not very 



close to T c . 
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Now let us put q = and expand in q . Since the constant term is canceled out by other 
diagrams, one has to look only at terms of order q 2 . One finds 



Z(q ,O) = -8X 2 v 2 q 2 J 



dp 



f(u p ) + f(n p ) + i , f(u P )-f(n p ) 



(Q p + Up) 3 (tt p - uj f 



(5.12) 



The first term in the square brackets gives rise to an integral that is dominated by the 
infrared, i.e., by p ~ m a , and so is completely analogous to the coefficient of p 2 above. The 
coefficient of q 2 coming from the first term is hence of order 0(X 1 ^ 2 ). In contrast, the integral 
of the second term in the square brackets is dominated by p ~ T. For such p one can write, 
approximately, 

n P -u p =^, / ( Wp )_ / (n p ) = _^^ l (5.i3) 



which gives 



8AV 2 f dp df(p) 4AV 2 



Substituting m 2 = 2At> 2 , we finally find 

S(g ,g) = -^ 2 - (5-15) 

Since T and f are both of order t> , the coefficient in front of q 2 is of order O(X ), in contrast 
to that of q 2 . The velocity of pions in our model is 

m2= ( i+ S) ' T>m - ( 5 - ie ) 

It is instructive to compare this result to the one obtained in the framework of the chiral 
perturbation theory: u — 1 — 0(T A /f£) |TJ| |2D|. One should bear in mind that the regime 
in which the chiral perturbation theory result applies corresponds to T m a . In our linear 
sigma model we consider a different regime: T ~ v ^> m a ~ y/Xv. In the chiral perturba- 
tion theory T/f n serves as an expansion parameter. In our weak-coupling calculation the 
expansion is in A while the T dependence is included to all orders in T/v. 

Now we need to check that this coincides with f 2 /X I5 . Recall that our definition of f 2 is 
as follows: if we define (p a = ir a /v, then 

dx e- iq - x ( V a (t, x)<p b (t, 0)) = — . (5.17) 



p q 2 

The left hand side is proportional to the ir Matsubara propagator, summed over q , 

( 5 - 18 ) 

go u 

When \q\ the dominant term in the sum is the one with go — 0, therefore, 

f 2 = v 2 . (5.19) 
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Taking into account one-loop graphs, as we have seen, will change f 2 only by an amount of 
order 0(A 1/2 ). 

Now to compute the susceptibility X I5 we need to turn on a chemical potential coupled to 
a broken charge. There are N — l broken generators. Let us consider the one that transforms 
a and tt\ into each other. The change of the Lagrangian when the corresponding chemical 
potential is turned on is 

u 2 

5C = M^kV - (rdoTn) + y [(v + a) 2 + u 2 } . (5.20) 

The susceptibility can be computed in Matsubara formalism. There are two contributions: 
one from the fi 2 term in 5C, the other from the bubble graph: 



2 1 f dp . . _ 2 p n T 2 

= v H — — / - — — sinn — — = v H . 

2T J (2vr) 3 2T 3 



2 

"2^2 



(5.21) 



The interpretation of this formula is rather direct: v 2 is the contribution of the condensate, 
and T 2 /3 is the contribution from the free gas of a and tti- The square of the velocity of 
pions, Eq. (|5.16|) , is equal to the ratio of f 2 in Eq. (|5.19|) and X I5 in Eq. ( |5.21|) , which is 



what we need to verify. We did not, however, attempted to turn on an explicit symmetry 
breaking and verify, e.g., the GOR relation at finite temperature. Such a calculation should 
be straightforward. 



VI. CONCLUSION 

Our goal has been to demonstrate that the dispersion relation of pions can be expressed 
in terms of quantities obtainable from equal-time or static correlation functions. The precise 
relation is given in Sec. [TI|. Our result enables one to find the real part of the pion dispersion 
relation on the lattice. However, it does not enable one to compute the imaginary part, 
which characterizes the damping of pion modes. 

Nowhere in our treatment did we assume any condition on the temperature (except T < 
T c ), as long as we stay infinitesimally close to the chiral limit. Thus, if m q is infinitesimally 
small, our result applies to all temperatures smaller than the temperature of the chiral 
phase transition T c . However, at any fixed (small) value of m q , our results do not apply 
at temperatures too close to critical, where the pion screening mass becomes of the same 
order as the screening mass of the sigma meson. The width of this region near T c shrinks to 
zero as a power of m q (more precisely, as m^ 135 ). Our treatment must fail there because, as 



explained in Sec. IV, the sigma boson also needs to be included into the hydrodynamic theory. 



However, the scaling of different quantities in this temperature region can be determined 
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using scaling and universality arguments, as discussed in Ref . . In the Appendix of the 
present paper we derive some additional interesting scaling properties omitted in Ref. M. 

Finally, assuming m q to be very small, as T — > T c , the chiral condensate {ipip) — > 0. As 
shown in Ref. |J, this implies that also / — > (although with a slightly different critical 
exponent). On the other hand, the axial isospin susceptibility X I5 becomes degenerate with 
the (vector) isospin susceptibility at T c , where both remains finite. One concludes from Eq. 
( p. 2D that the pion velocity u tends to zero as one approaches the critical temperature. In 
fact, it can be shown that u approaches zero faster than the divergence of the screening mass 
m, so the pion pole mass m p = um goes to zero as T — > T c ||. 
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APPENDIX A: DYNAMICAL SCALING AT THE CHIRAL TRANSITION 

In this appendix we provide a simple derivation of the dynamical critical exponent z 
characterizing critical slowing down near the QCD chiral phase transition. It is based on a 
scaling argument similar to the one used in Ref. [18| to find z in antiferromagnets. 

For simplicity, here we consider the chiral limit m q = 0. Scaling and universality argu- 
ments presented in Ref. |§ predict that the velocity u vanishes as T approaches T c from 
below. The quantity u 2 scales as u 2 ~ 

t (d-2)u^ where t _ _ t)/T c . Since the inverse 
correlation length of the order parameter ipip, i.e., the static screening mass of the sigma 
particle m a , scales as m a ~ t v , we conclude that for tcl 

u 2 ~ m d ~ 2 . (Al) 

The fact that u —>■ at T c means that the dispersion relation ceases to be linear. Moreover, 
the effect of damping also becomes important; in other words, we expect the real and 
imaginary parts of u to become comparable. In such a situation on cannot refer to u as a 
quasiparticle energy. Rather, it is a characteristic frequency, or inverse relaxation time, of 
a mode of a given wave number p. Scaling hypothesis dictates that the relation between 
p and u> should be homogeneous, i.e., u ~ \p\ z . The dynamical scaling exponent z is, in 
a generic case, new exponent independent of the static exponents, e.g., v and 77. However, 
in the case of QCD (similar to the case of an antiferromagnet |6|, [L8|), since the dispersion 
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relation of pions is given in terms of static quantities only, it turns out that z, as one would 
expect, can be derived from static scaling laws only. 

To determine z we observe that the dispersion relation uj ~ \p\ z applies at scales m a <C 
\p\ <C T. At softer scales, |p| <C m a , the dispersion relation is still linear, uj = u\p\. 
Requiring that the two expressions for uj match at |p| ~ m a we find 

m^~Mm (J ~m^ /2 , (A2) 

where to obtain the last scaling relation we used Eq. ( |Al|) . Therefore, at the QCD chiral 
phase transition the dynamical critical exponent is given by 

- = ~ (A3) 

(see also an alternative derivation in Refs. || |18|). 6 

A similar argument can be applied to determine the scaling of the diffusion coefficient 
D' — D + K2, characterizing the pion damping in the chiral limit [Eq. (|4.22b )j. The small 



momentum expansion of uj is u\p\ — ^D'p + . . . . The scaling hypothesis dictates that all 
terms in this expansion become of the same order in magnitude when \p\ ~ m a (otherwise, 
another scale, e.g., u/D', appears in addition to the inverse correlation length m a ). This 
requires 

D' ~ mi d - 4 ^ 2 ~ t^ d -^l\ (A4) 
Thus D' diverges near the chiral phase transition in QCD (d = 3). 
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